The Influence of Magnetic Anisotropy on the Kondo Eff'ect and Spin-Polarized 
Transport through Magnetic Molecules, Adatoms and Quantum Dots 
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Transport properties in the Kondo regime of a nanosystem displaying uniaxial magnetic anisotropy 
(such as a magnetic molecule, magnetic adatom or quantum dot coupled to a localized magnetic 
moment) are analyzed theoretically. In particular, the influence of spin-polarized transport through 
a local orbital of the system and exchange coupling of conduction electrons to the system's magnetic 
core on the Kondo effect is discussed. The numerical renormalization group method is applied to 
calculate the spectral functions and linear conductance in the case of the parallel and antiparallel 
configurations of the electrodes' magnetic moments. It is shown that both the magnetic anisotropy 
Eis well as the exchange coupling between electrons tunneling through the conducting orbital and 
magnetic core play an important role in formation of the Kondo resonance, leading generally to 
its suppression. Specific transport properties of such system appear also as a nontrivial behavior 
of tunnel magnetoresistance. It is also shown that the Kondo effect can be restored by an external 
magnetic field in both the parallel and antiparallel magnetic configurations. 



I. INTRODUCTION 

Growing trend towards building ever more efficient 
and smaller electronic devices inevitably draws the re- 
searchers' attention to nanoscopic hybrid systems. In this 
respect, single atoms or molecules seem to be promis- 
ing prospects, as their incorporation into electronic nan- 
odevices allows for developing novel systems capable of 
performing strictly imprinted functions^"— among which 
information storage is of key interest i2r— Consequently, 
due to recent advances in experimental techniques en- 
abling to address transport through individual atoms and 
molecules, both natural as well as artificial (quantum 
dots) systems exhibiting magnetic anisotropy, such as 
magnetic atoms of spin S > 1/2 (i.e. Fe, Co or Mn)^iii"— 
or single-molecule magnets (SMMs)fi^"— have become 
the object of intensive studies. 

It has been suggested, and in the case of magnetic 
adatoms also experimentally proven;^ that magnetic 
state of such systems can by controlled by the use of 
spin-polarized currents^"— or spin bias;^ This practi- 
cally means that the system's magnetic moment can be 
switched between two metastable states of minimal en- 
ergy by only applying an electric/spin current pulse of 
a proper amplitudcj^ Furthermore, if attached to two 
metallic nomnagnetic electrodes a SMM can act as a 
spin filter<^"— If, however, coupled to electrodes char- 
acterized by unequal spin polarizations, the molecule 
can reveal transport characteristics typical of a spin 
diodci^ Most of these results have been obtained in the 
limit of weak coupling between a SMM and reservoirs of 
spin-polarized electrons. Nevertheless, in some situations, 
when mixing of localized electron states responsible for 
transport properties of the molecule and extended elec- 
tron states in electrodes is significant, such an assumption 
not necessarily has to be correct. 



In the strong coupling regime the electronic correla- 
tions can lead to an additional resonance in the den- 
sity of states at the Fermi level of electrodes, known 



as the Kondo- Ahrikosov-Suhl 



resonance 
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end of the 1990s, the presence of the Kondo effect 
has been successively demonstrated in a large vari- 
ety of nanoscopic objects like quantum dotsj^^"— mag- 
netic adatoms nanotubes^ and different types of 
molecules: Co(II)-based coordination complexesj^ di- 
vanadium moleculeS)^ and Cgo molecules attached to 
gold^ or ferromagnetic nickel electrodes4i However, in 
the case of nanosystems characterized by large spins the 
prominent role of the magnetic anisotropy in formation 
of the Kondo effect has been experimentally established 
only very recentlyjiii^ It turned out that the Kondo 
effect can be then tuned by changing both the orien- 
tation (e.g. by controlling the adatom's local environ- 
ment^i^) and magnitude (e.g. by mechanical straining of 
the molecule^) of the magnetic anisotropy. Moreover, 
Parks et al^ were able to tune the anisotropy constant 
continuously and to modify accordingly the energy spec- 
trum underlying the Kondo state. As a result, they man- 
aged to observe a crossover from the fully screened to 
underscreened Kondo effect. It is worthy of note that 
more recently electric field control of magnetic anisotropy 
has been experimentally established for a SMM embed- 
ded into a planar three-terminal device^S Although a 
few theoretical works focused on transport related issues 
in SMMs in the Kondo regime have been already pub- 
lishedji^"— experimental evidence of the Kondo effect in 
transport through SMMs has been found only very re- 
cently.i^ 

The earlier works on the Kondo phenomenon in trans- 
port through SMMs have been primarily focused on the 
role of transversal magnetic anisotropy, and hence also on 
the role of quantum tuimeling of the SMM's spin in the 



formation of the Kondo state. It has been shown that 
the interplay of quantum tunnehng and spin exchange 
processes between the molecule and tunneling electrons 
may result in the pseudo-spin 1/2 Kondo efFect4^i^ Fur- 
thermore, it was soon realized that when even a mod- 
erate transverse magnetic field is applied, any qualita- 
tive difference between the mechanisms of the Kondo ef- 
fect for molecules with half- and full-integer spins cease 
to exist4^ Mapping of the Anderson-type Hamiltonian 
describing a SMM onto the spin- 1/2 anisotropic Kondo 
Hamiltonian^ led to the conclusion that, depending on 
whether the molecule's total spin is reduced or aug- 
mented in the singly charged state, the coupling in the 
Kondo Hamiltonian is antiferromagnetic or ferromag- 
netic, respectively. In the former case the Kondo effect 
is revealed, whereas in the latter one no resonance at 
the Fermi level is present due to renormalization of the 
transverse coupling to zero. In addition, the Kondo effect 
is expected to oscillate as a function of the magnitude of 
transverse field due to the Berry-phase periodical mod- 
ulation of the tunnel splitting42. Finally, the nonequilib- 
rium spin dynamics of a SMM, triggered by a sudden 
change in the magnetic field amplitude has been stud- 
ied, with the main emphasis on the time evolution of the 
Kondo screening 1^ 



Since physical mechanisms governing Kondo correla- 
tions in spin-polarized transport through nanoscopic sys- 
tems exhibiting the magnetic anisotropy, such as mag- 
netic adatoms or SMMs, are still at the early stage of 
research, the objective of the present paper is to provide 
further insight into the problem. In particular, we investi- 
gate how magnetic anisotropy affects the system's trans- 
port characteristics such as conductance and tunnel mag- 
netoresistance (TMR) in the linear response regime. To 
properly describe the transport properties in the strong 
coupling regime, we employ the Wilson's numerical renor- 
malization group (NRG) approachj^2i^i^ This method is 
known as very powerful and essentially exact in solving 
quantum impurity problems. 



The paper is organized as follows. In Sec. II we describe 
the model Hamiltonian used in calculations and provide a 
brief introduction to the NRG method. Numerical results 
and their discussion are given in Sec. HI, where we ana- 
lyze the spectral functions of the orbital level as well as 
the conductance and tunnel magnetoresistance (TMR) in 
the linear response regime. The above quantities are an- 
alyzed as functions of the orbital level position, strength 
and type of exchange coupling, and the anisotropy con- 
stant. In addition, we also discuss the effect of external 
magnetic field. Finally, the summary and conclusions are 
given in Sec. IV. 
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Figure 1. (Color online) (a) Schematic representation of the 
system under consideration. The system consist of two ferro- 
magnetic electrodes to which a magnetic quantum dot (MQD) 
exhibiting magnetic anisotropy is attached. As the MQD one 
can conceive either a magnetic adatom (i.e. Fe, Co, Mn) (b), 
semiconductor quantum dot coupled to a magnetic moment 
(c), or a single- molecule magnet (SMM) (d) - here, only the 
magnetic core of the Fe4 molecule^ is schematically depicted. 

II. THEORETICAL DESCRIPTION 
A. Model 

We consider a generic theoretical model that allows 
for capturing essential features of quantum objects such 
as magnetic adatoms, quantum dots coupled to localized 
magnetic moments, and SMMs, see Fig. [T] It is assumed 
that electronic transport takes place via a single local 
orbital level (OL) of the system (conducting orbital of a 
SMM, adatom or quantum dot), which is coupled to elec- 
trodes. Moreover, the OL is also exchange-coupled to the 
corresponding magnetic core. Without loss of generality, 
we will henceforth refer to the systems under investiga- 
tion as magnetic quantum dots (MQDs). 

The total Hamiltonian of a MQD coupled to external 
leads can be written as 

Ti = "HmQD + Hlcads + Htun- (1) 

The first term represents the MQD and has the fornt^i"— 
Hmqd = - DSl + ^ en„ + U n^n\^ 

-Js-S + B,{S,+s,), (2) 

where D stands for the uniaxial anisotropy constant of 
the MQD, while Sz denotes the zth component of the 
MQD's internal spin operator S. Since in the present pa- 
per we focus only on systems displaying magnetic bista- 
bility. the anisotropy constant is assumed to be positive 
{D > 0). Furthermore, ricr = cj.c^ is the OL occupa- 
tion operator, where cJ(Cg.) creates (annihilates) a spin- 
(T electron of energy e in the OL. The Coulomb energy 
of two electrons of opposite spins occupying the OL is 
given by U . The penultimate term of Eq. ^ accounts 
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for exchange couphng between the magnetic core of a 
MQD and the spin of an electron in the OL, represented 
by s = \ Y.aa' cl(T„„,c^,, whcrc (T EE (cr^ , CT?' , (T^ ) is the 
Pauh spin operator. The J-couphng can be either of fer- 
romagnetic (J > 0) or antiferromagnetic (J < 0) type. 
FinaUy, the last term of Eq. ([2]) describes the Zeeman 
interaction of the MQD with an external magnetic field 
B — (0,0, Bz) oriented along the easy axis of a MQD. 
Note that we put here gfiB = 1- 

The ferromagnetic metallic electrodes, to which a 
MQD is coupled through the OL, are characterized by 
noninteracting itinerant electrons with the dispersion re- 
lation e^^, where q indicates either left {q = L) or right 
{q = R) electrode, k denotes a wave vector and cr is a 
spin index of an electron. Thus, the leads' Hamiltonian 
is given by 



Hi 



cads 



E EE 

q=L,R k a=t,l 



(3) 



with a^lj. (a^^) being the relevant electron creation (an- 
nihilation) operator. At this point, it should also be men- 
tioned that in the present paper we limit the discussion 
to coUinear (parallel and antiparallel) configurations of 
electrodes' magnetic moments. Furthermore, the MQD's 
easy axis is assumed to be collinear with these moments 
as well. Finally, electron tunneling processes between the 
MQD and electrodes are included in the term Htun, 



-Ht 



E E E ^fco-^fc^ 

q=L,R k cr=-|',i 



H.c. 



(4) 



where T? denotes the tunnel matrix element between 
the OL and the gth lead. 

In the linear response regime, it is numerically con- 
venient to introduce the following canonical transforma- 
tion 
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''ka \ _ 
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^k(j 



^k(7 



fca ^ka 



'■ka 



(5) 



where, Vka = ^\Tt? + \Tt?, and the label e (o) de- 
notes the even (odd) combination of the leads operators. 
Such a rotation in the space of the left-right electron 
operators results in separation of the total Hamiltonian, 
Eq. Ill, into two independent parts. The first one in- 
volves the OL coupled to a single electron reservoir de- 
scribed by the even linear combination of the leads' elec- 
tron operators, a|^, while the other one is related with 
non-interacting electrons decoupled from MQD and de- 
scribed by the odd operators a,%„. Consequently, the tun- 
neling Hamiltonian can be written as 



Htun -EE 

k a=\.\. 



Vko 



(^ka^^a 



■ cia 



ka 



(6) 



with Vka being effective OL-lead tunneling matrix ele- 
ments. In the following we assume that the full spin- 
dependence is taken into account through the matrix el- 



ements Vkc 
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For simplicity, we assume a flat con- 



duction band in the interval [— so that Pcr(a;) 



p = '^ith W representing the cut-off energy of the 
system and W = 1 taken as the energy unit. Finally, 
the energy dependence of Vka is neglected, Vka ^ Va .— 
Under these circumstances, the overall effect of the ferro- 
magnetic reservoir on the MQD is completely determined 
by the hybridization function Ta, 



^P\Va\ 



(7) 



B. Method of calculations 



In order to determine transport properties in the 
strong coupling regime, we use the Wilson's numerical 
renormalization group methodi^i^i^ The NRG tech- 
nique consists of logarithmic discretization of the conduc- 
tion band (with a discretization parameter A > 1) into 
intervals [A-("+i)H^, A-"H^] and [-A-"VF, -A-("+i)VF] 
for ?7 = 0, 1, 2, 3, . . ., which allows for resolving transport 
properties on energy scales logarithmically approaching 
the Fermi level. After having discretized the conduction 
band, such a model is mapped onto a semi-infinite chain, 
whose first site is coupled to the impurity (in our case 
the MQD). The Hamiltonian then reads^i^ 

n = Hmqu + V ^ [^'^■^o<T + fLca] 

oo 

+ E/ E/ ^"^[fna-fn+la+fn+lafna]- (8) 
n=0 (T=t,i 

The operators {f^^) correspond to the nth site 
of the Wilson chain, with exponentially decaying hop- 
ping matrix elements t„ between neighboring sites of the 
chain<^ As a consequence, by adding consecutive sites, 
one is able to access transport at lower and lower en- 
ergy scales. In this way the method generally provides a 
non-perturbative description of the crossover from a free 
magnetic impurity at high temperatures to a screened 
spin at low temperatures.— The Hamiltonian ([8]) can 
be solved iteratively by adding consecutive sites of the 
chain. This procedure allows for resolving static and dy- 
namic properties of the system at energy scale A~"/^, 
with n being a given iteration. 

Since the NRG calculations may in general pose a se- 
rious numerical challenge, it becomes essential to take 
advantage of as many available symmetries of the sys- 
tem's Hamiltonian as possible. To efficiently address 
the present problem, we have employed the fiexible 
density-matrix numerical renormalization group (DM- 
NRG) codc)^ which can exploit an arbitrary number of 
both Abelian and non-Abelian symmetriesi^ In the case 
under discussion, the C/charge(l) x t4pin(l) symmetry of 
the model has been used, so that the zth component of 
the total spin. 
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n=0 



(9) 
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where Si = Sz + Sz , as well as the total charge 



(7 n—0 \ (7 



(10) 



served as quantum numbers according to which the states 
of the Hamihonian were classified during computation. 
Finally, the discretization parameter A = 1.8 has been 
taken in calculations, and we have kept 2000 states after 
each step of the iteration. 

The central quantity of interest is the OL spin- 
dependent spectral function^'^ 



1 



Im((cJ4)) 



(11) 



where ((Co.|cJ.))L denotes the Fourier transform 
of the retarded Green's function ((Co.|cJ.))( = 
— z6'(i)({Cg.(t), c^. (0)}) of the orbital level. Some technical 
details concerning calculation of the spectral function 
can be found in Appendix 1X1 

Having found the spectral function, one can determine 
the linear response conductance g from the Landauer- 
Wingreen- Meir formula^"^ which at T = yields 



A^uj = 0) (in units of ^), (12) 



with r^(^) = ^{1 ± P) and r^(j^) = ^(1 ± P) for 
the parallel magnetic configuration of electrodes, while 



5(1 ± P) and r^^^j = E(l ip P) for the antipar- 



allel one. We assumed that both electrodes are made 
of the same material and P denotes their spin polar- 
ization. The effective coupling between the MQD and 
the reservoir for the parallel magnetic configuration is 
^^{i) ~ -'^(1-'-^)' '^hile for the antiparallel one Fi^^^ = F, 



with F = (F| + F;)/2. Note that in the case of left-right 
symmetric systems, is independent of a in the an- 

tiparallel configuration, and hence the system behaves 
effectively as coupled to nonmagnetic leads. As a result, 
the linear spin-resolved conductance (in units of ^) in 
the two magnetic configurations is given by 



3m) = o(i±^)r4(i)(^ = o)' 



9n\)-^i^-P')^Afa)i^^o), 



(13) 



with Aa^^^ being the spectral function in the respective 
magnetic configuration. 

III. NUMERICAL RESULTS AND DISCUSSION 

Transport characteristics of the system, such as the 
linear-response conductance g and tunnel magnetoresis- 
tance (TMR) , have been numerically calculated for a hy- 
pothetical MQD characterized by the spin number S = 2. 



Since the key feature of the transport regime under dis- 
cussion is the presence of the Kondo resonance, it is con- 
venient to introduce the Kondo temperature Tk as the 
most relevant energy scale for the system, to which other 
energy parameters will be related, if necessary. In the 
present work, the Kondo temperature Tk is estimated 
from the half-width at half-maximum of the Kondo res- 
onance in spectral function at T = for J = and 
P ~ 0^1^ We note here that, for the sake of simplicity, 
we assume fee = 1, i.e. temperatures are also given in 
units of energy. As a result, for the parameres used in 
Fig.H we get Tk/W w 5 • lO"'^ (Tk/F « 0.022). 

Before presenting and discussing numerical results, it 
is worth recalling that the Kondo effect appears as a re- 
sult of spin exchange processes in the OL due to its strong 
coupling to electrodes. This coupling is described here by 
the effective hybridization parameter F, which introduces 
the Kondo temperature Tk as the relevant energy scale. 
However, the model considered provides also another in- 
dependent physical mechanism (channel) through which 
the Kondo effect can be modified, i.e. the exchange inter- 
action J between an electron in the OL and the MQD's 
magnetic core. Therefore, one can expect the ratio of the 
J-coupling and the Kondo temperature Tk to be the key 
parameter controlling whether the Kondo resonance will 
appear or not. Indeed, we show below that the electronic 
correlations between the OL and electrodes effectively re- 
sult in formation of the Kondo resonance only if | J| < Tk- 
In order to see how the parameters of the system, es- 
pecially the J-coupling and the magnetic anisotropy D, 
influence the transport properties, we first calculate and 
discuss the relevant spectral functions. 



A. Spectral functions 

In the following we will analyze dependence of the spec- 
tral functions of the OL on some essential parameters of 
the system. Figure [5] shows the normalized spin-resolved 
spectral functions T:Ta-A^{Lu) as a function of the OL en- 
ergy e in the antiparallel and parallel magnetic configura- 
tions for I J|/Tk ~ 2. In the antiparallel configuration the 
spectral functions for both spin components are equal. In 
the parallel configuration, on the other hand, the main 
contribution comes from the spin-up electrons which are 
the majority ones. In both configurations and for both 
types of exchange coupling, the spectral functions show 
clear resonances associated with degeneracy of the neigh- 
boring charge states - compare the boundaries between 
regions corresponding to different Q, where Q denotes 
the average number of electrons occupying the OL. In 
the singly-occupied regime, Q = I, the Kondo effect due 
to hybridization of the OL spin with the conduction elec- 
trons of the leads should be present. However, there are 
two ingredients that may generally suppress the Kondo 
resonance: the exchange coupling J and the exchange 
field due to the presence of ferromagnetic leads. Since the 
results in Fig. [2] are shown for | J|/Tk ~ 2, only the remi- 
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Figure 2. (Color online) Normalized spin-resolved orbital level (OL) spectral functions, -kV crAdiS), shown as a function of 
the OL energy e in the antiparallel (a)-(b) and parallel (c)-(f) magnetic configurations. Top panel corresponds to the case of 
ferromagnetic ( J > 0) coupling between electrons in the OL and MQD's magnetic core, while the bottom one presents results 
for the antiferromagnetic coupling. The variable Q denotes the average number of electrons occupying the OL. The parameters 
are: U/W = 0.3, D/U 1.7 ■ 10-* (D/Tk = 0.1), T/U 0.075, | J|/r ^ 0.044 (| JI/Tk « 2), B, = 0, and P = 0.5. Note that 
the spectral functions are presented in a logarithmic scale. 



niscent of the Kondo effect can be observed. The suppres- 
sion of the Kondo resonance is especially visible for the 
antiferromagnetic coupling (J < 0). Moreover, while in 
the case of ferromagnetic coupling ( J > 0) some residual 
Kondo resonance can still be visible for the energy corre- 
sponding to the particle-hole symmetry (e = — C//2), the 
resonance is practically absent for the antiferromagnetic 
coupling, though some side resonances appear. The ori- 
gin of these additional features will be discussed further 
in the text. 

Let us now analyze how the shape of spectral functions 
in the Kondo regime evolves when the exchange interac- 
tion between electrons in the OL and the MQD's core is 
turned on gradually, see Figs. [3][a)-(f). Note that when- 
ever we consider behavior of the system in the particle- 
hole symmetric point, i.e. for e = —U/2, only the range 
of positive energies is presented. As one might expect, 
the behavior of the system for small values of |J|, i.e. 
1^1 ^ Tk, resembles that of a single-level quantum dot, 
and a well defined and pronounced Kondo peak in the 
antiparallel configuration of the electrodes' magnetic mo- 
ments is observed for oj < Tk, see Figs. [2I^a)-(b). In the 
parallel configuration, on the other hand, spin-dependent 
coupling to the electrodes acts as an effective exchange 
field^S leading to spin splitting of the OL. This in turn 
results in suppression of the Kondo effect, except for the 
particle-hole symmetric point, e = —U/2, as shown in 
Figs.[3]Jc)-(f); where the effective exchange field vanishes. 



In the antiparallel configuration, the resultant coupling 
is the same for the spin up and spin down when the sys- 
tem is left-right symmetric. Consequently, there is no ex- 
change field and we observe a well-pronounced peak in 
the spectral function at the Fermi level also outside the 
particle-hole symmetric pointi^Sr— 

The height of the Kondo peak becomes reduced with 
the increase of |J|, and for |J| ^ Tk the peak almost 
completely vanishes. We note that the disappearance of 
the resonance is faster in the case of the antiferromag- 
netic coupling (J < 0), see Figs. [3{b,d;f). Furthermore, 
as the J-coupling grows, some additional features in the 
spectral functions emerge. Apart from the Hubbard peak 
originating from the Coulomb repulsion of two electrons 
in the OL level, there are two additional resonances for 
J > 0, marked as dashed lines A and B in Fig. [3][a), and 
one resonance for J < 0, line C in Fig. [2tb). Interest- 
ingly, position of one of the two resonances for J > 
remains roughly independent of energy (line A) , whereas 
the other resonance moves towards larger energies as J 
increases (line B), see Fig.[31[a). 

Some insight into physical origin of the resonances A, 
B and C can be gained by considering the lowest energy 
states of a free-standing MQD with one extra electron in 
the OL, see Figs. [3l^g,h).— First, we note that the con- 
sequence of exchange interaction between an electron in 
OL and magnetic core is a decomposition of the molecular 
magnetic states into two spin-multiplets, corresponding 
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Figure 3. (Color online) Normalized spin-resolved orbital level 
(OL) spectral functions shown for different values of the 
exchange parameter J in the case of (a,c,e) ferromagnetic 
(J > 0) and (b,d,f) antiferromagnetic (J < 0) coupling, and 
for e — —U 12. The spectral function displays two additional 
resonances for J > 0, marked by dashed lines A and B, and 
a single resonance for J < 0, marked by a dashed line C. 
The bottom panel (g,h) illustrates the dependence of several 
lowest energy states of a singly occupied MQD on the pa- 
rameter J (thick lines), and the corresponding probabilities 
of finding the electron in a certain spin-state for j± (thin 



lines): o? = {A+f = (1^)2, while = {M+f = (A-)^ [for 
details see Eq. (|14|l and the paragraph below it). Remaining 
parameters as in Fig. [51 



to 5 + 1/2 and 5—1/2. In addition, the sign of the cou- 
pling parameter J determines which of the two multiplets 
has lower energy. Since we focus exclusively on the case 
of T = 0, it is justified to take into account only the rele- 
vant low energy states in both spin-multiplets. These are 
presented in Figs. [3]Jg,h), where the zero- field {Bz — 0) 
energy of states |± |), and |± is presented as a func- 



tion of the couphng parameter J. Tire superscript ± at 
the states | ± |)* is used to distinguish between states 
of higher (-I-) and lower (— ) energy. We note that the 
state |± I) belongs to the spin-multiplet corresponding 



to 5 + 1/2, while for J > the state |±|) 



belongs to the 



multiplet S 
5—1/2 (note we assumed 5 



1/2 and the state |± |)+ to the multiplet 



is opposite, i.e. |± |) belongs to the multiplet 5 — 1/2 
whereas the state |± I)"*" to the multiplet 5 + 1/2. In the 
absence of external magnetic field the MQD's states of 
interest take the form 



Si 
Si 

Si 



±f) = lta))oL®|±2)corc, 

,i;)oL® 1-1) 

)oL 



-3/2 I 



8^3/2lt)c 



core 1 



(14) 



= +f) = A^3/2lt)0L® |+l)core 
+ B:};3/2|;)oL ® |+2)core, 

with |»)oL(corc) denoting the spin state of OL (mag- 
netic core). The coefficients A^^ = cxp[i(f>^] and 
= exp[i7^] can be regarded as effective Clebsch- 
Gordan coefficients, which are nontrivial functions of the 
system's parameters J and Here, (j)^^ and 7^ are rel- 
evant phase factors. As above, the superscript ± is used 
to distinguish between states of higher (+) and lower 
(— ) energy. It is worth emphasizing that, to some ex- 
tent, the situation under consideration is similar to the 
case of a quantum dot subjected to an external magnetic 
field which leads to splitting of the Kondo resonance . 
However, in the case of a simple quantum dot in a mag- 
netic field there are two energy levels, whereas the energy 
structure of a MQD is much more complex, even in zero 
field. Moreover, except for the states |± |), all MQD's 
states for Q = 1 correspond to an electron in the OL be- 
ing in the superposition of spin-up and spin-down states. 

Let us now discuss cotunneling processes leading to 
the resonances for ferromagnetic (FM) exchange coupling 
(J > 0), the top panel of Fig. [31 We emphasize that now 
we consider the regime of large J (i.e. J > Tk and J < T), 
where the zero-energy resonance is absent and only side 
resonances (lines A and B) appear. Assume that initially 
the molecule occupies the state | — |). Due to spin-flip 
cotunneling processes, the MQD can be excited to one of 
the two states: |— |)~ (resonance A) and |— I)"*" (reso- 
nance B), see Fig. [Sjg). Taking into account the energy 
spectrumiS^ one can estimate the corresponding energy 
gaps for I J| ;2> -D as (exact formulas in Appendix [B|) 



FM 



25i:» 



1 - 



2i\J\-D) 



i2S+mj\-2D) 



A 



FM 



25 + 1 



(15) 



Jl 



2). For J < 0. the situation 



Thus, the resonance corresponding to the line A in 
Fig. [Sjja) is related to transitions characterized by the 
energy gap A™, while the resonance indicated by the 
line B is associated with the gap Af^. Moreover, one 
can note that Af ^ depends linearly on J, whereas A™ 
only slightly changes with J. 

The picture presented above for J > alters only 
slightly when the exchange coupling changes to the an- 
tiferromagnetic (AFM) one (J < 0), see Figs. [3Ub.d,f). 
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Figure 4. (Color online) Similar as in Fig. |3] but now the 
dependence of the normalized spin-resolved spectral func- 
tion on the uniaxial anisotropy constant D is presented for 
\J\/Tk ~ 2. a monochromatic color scheme is used here to 
highlight how the position of the resonances changes with 
varying D. 



Position of the spin-multiplets 5 -f 1 /2 and 5 — 1 /2 is now 
interchanged when compared to that for J > 0, Fig.[3^h). 
When the MQD is initially in the state |— |)^, then the 
spin-flip cotunneling processes can excite the MQD to one 
of the states: |-|) and |-|) + . Moreover, for | J| > D, the 
energy gaps associated with these transitions are roughly 
equal, 



AFM 



A 



AFM 



2S+1 



J|. 



(16) 



As a consequence, only the resonance denoted by the line 
C is then visible in Fig.[31[b). Physical origin of the reso- 
nances in the parallel magnetic configuration, Figs, ^c)- 
(e) and (d)-(f), can be accounted for in a qualitatively 
similar way. Finally, weak vertical lines visible in Fig. [3] 
for oj = U , correspond to the resonance between singly 
occupied and empty OL. 

As shown above, the exchange coupling of electrons in 
the OL and magnetic core modifies the energy spectrum 
of a MQD, and hence affects the Kondo effect. However, 
it has been demonstrated experimentall y^^i^°i^^ that the 



energy spectrum can also be modified by changing the 
anisotropy constant D. Variation of the OL spectral func- 
tions with D is shown in Fig. U] for a specific value of 
|>/|/7k ~ 2. Since \ J\ > Tk, the zero-energy Kondo reso- 
nance is practically absent, and only side resonances are 
visible. These resonances - marked in Figs. |31[a)-(b) by 
letters A, B and C - correspond to the relevant reso- 
nances in Figs. [31Ja)-(b). 

Especially interesting seems to be the case of J > 0, 
where two resonances (A and B) emerge at Z? w Tk. By 
a closer inspection of Fig. IH^g) one finds that this takes 
place when the condition D ^ J/2 is satisfied. However, 
it should be noted that according to our definition of the 
Kondo temperature, the relation D = J/2 Ri Tk is only 
coincidental, and thus valid just for the current set of 
parameters. Moreover, in the limit D 3> |J| the energy 
gaps introduced above can be estimated as 



FM 



(25 + l)V' 
16(25'- 1)1?' 

(25 + l)V2 



SJ- 



A FM 
^2 



(25- i)i:» 



(17) 



16(25- I)!?' 



so that A\ is nearly constant with respect to D, 
whereas A™ depends almost linearly on D. At first 
glance, the latter result appears to contradict the shape 
of the resonance B seen in Fig.lDJa), so that consideration 
of the energy spectrum becomes apparently insufficient. 
In order to account for this disparity, we must also take 
into account the explicit form of the MQD's states con- 
tributing to the resonances, see Eq. ([14]). For D > \ J\/2, 



we have 



-y > {A-y and (A+)^ > 



and as the 



at some point we get 
(B+)2 « 0, Fig. Hg), 

- |t)oL®h2)coro and 



magnetic anisotropy D increases 
(B-y = {A+y « 1 and (A-)^ = 
which basically means that |— 
|— I)"*" ~ |i)oL® |- l)core- Li Other words, when D > \ J\, 
the cotunneling processes associated with transitions be- 
tween the states |— |) and |— |)^ are allowed, while those 
between the states | — f ) and | — f)"*" are forbidden. Con- 
sequently, the resonance marked with a vertical dashed 
line in Fig. IHb) actually represents two independent res- 
onances: B for D <\J\ and A for Z) > | J|. 

To complete the discussion of spectral functions, we 
note that the presence and behavior of the resonance C 
for the antiferromagnetic J-coupling (J < 0), Fig. (H^b), 
can be explained in a way similar to that for J > 0, with 
the relevant energy gaps for D ^ \J\ estimated as 



AFM 



S\J\ 



A 



AFM 



(25- i)i:» + 



(25-H)V 
16(25-1)13' 

(25+l)V2 



(18) 



8(25- i)i:» 



However, only virtual spin-flip transitions between the 
states I — |)^ and | — |), and represented by the en- 
ergy gap A^^^ [see Eq. ([T5)) ]. are then possible. On the 
other hand, in the opposite limit, D ^ | J|, both types of 
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Figure 5. (Color online) Total linear conductance g = J^^. 5ct in the antiparallel (a,b) and parallel (c,d) magnetic configurations, 
and the spin-resolved linear conductance in the parallel configuration (e)-(h), presented as a function of the OL energy e 
for indicated values of the parameter J in the case of the ferromagnetic ( J > - top panel) and antiferromagnetic (J < - 
bottom panel) exchange interaction between electrons in the OL and magnetic core. All other parameters are as in Fig. [2] 



spin-flip cotunneling processes characteristic for the an- 
tiferromagnetic J-coupling (as discussed earlier) can in 
principle operate. Both these transitions are associated 
with the same energy gap given by Eq. ([T5|. which, un- 
like in the case of J > 0, results in formation of one 
resonance only. Fig. Hl^b) . 



B. Conductance in the linear response regime 

From the spectral function discussed above one can de- 
termine the spin-resolved as well as total linear conduc- 
tance, shown in Figs. [SJa)-(g). For \ J\ <^ Tk, the results 
well known for Kondo effect in a single-level quantum 
dot are recoveredi^Si^iZi In particular, for the antiparallel 
magnetic configuration an enhanced conductance occurs 
in the blockade regime (single electron in the OL, Q = 1), 
whereas in the parallel configuration only a sharp peak 
in the particle-hole symmetry point, e = —U/2, appears. 
More precisely, in the antiparallel configuration the lin- 
ear conductance (measured in the units of 2e^ /h) is given 
by g^^ = 1 — , while in the parallel configuration the 
conductance reaches unity for e = —U/2, = 1, with 
g,P/AP _ g^^^^ representing the total linear conduc- 
tancCf^iiZ^ Suppression of the Kondo effect for other val- 
ues of e in the latter case is a consequence of spin splitting 
of the OL due to an effective exchange field created by fer- 
romagnetic electrodes, as it has already been mentioned 
in the previous section. 

When strength of the J-coupling increases, the Kondo 
effect becomes gradually suppressed and the linear con- 
ductance in the blockade {Q = 1) region decreases 
In order to explain this dependence, we note that in 



the situation under discussion electrons are transmitted 
through OL that is exchange-coupled to magnetic core, 
see Eqs. (|T4|) . As a result, cotunneling processes respon- 
sible for the Kondo state are more complex than in the 
case of J = 0. First, the amplitude of the cotunneling 
processes becomes reduced as now the electron occupying 
the OL is in a superposition of the spin-up and spin-down 
states. Second, the J-coupling creates an energy gap be- 
tween the relevant states, which effectively suppresses the 
Kondo effect. 

Looking more carefully at the conductance curves in 
Fig. O one finds some difference between the cases with 
J > and J < 0. While for J < the Kondo peak van- 
ishes rapidly after | J| exceeds Tk, some remanent Kondo 
peak is still visible for J > 0. This disparity stems from 
different properties of quantum states taking part in the 
formation of the Kondo state for J > and J < 0. 
First, when | J| > Tk, the ground state energy of a singly- 
occupied MQD is lower for J < than for J > [compare 
Figs.[3]^g) and (h)]. Second, the ground state for J < is 
a superposition of spin-up and spin-down states, which 
is not the case when J > 0. As a result, the cotunneling 
processes driving the Kondo effect are more effective for 
J > than for J < 0; which leads to lower conductance 
for J < as compared to J > 0. 

From the previous subsection we know that the uni- 
axial magnetic anisotropy modifies electron states of a 
MQD, affecting thus the Kondo effect. This is shown ex- 
plicitly in Figs. [5] and [T] As one can note in Fig. [SI con- 
ductance in the parallel magnetic configuration is rather 
insensitive to the anisotropy constant D - regardless of 
the type of exchange coupling J, see the filled points 
in Fig. [7)Ic)-(d), and certain weak dependence on the 
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Figure 6. (Color online) Influence of the uniaxial magnetic 
anisotropy D on the total linear conductance g = for 
|J|/Tk ~ 2 and the exchange J-coupling of either (a,c) fer- 
romagnetic (J > 0) or (b,d) antiferromagnetic (J < 0) type. 
Other parameters are the same as in Fig. [21 

anisotropy constant D appears then in the particle-hole 
symmetry point e — ~U/2 (Kondo peak), see Figs.[71[a)- 
(b). This dependence, however, becomes weaker when 
the magnitude of exchange coupling increases. On the 
other hand, the role of anisotropy is more important in 
the case of antiparallel magnetic configuration (also for 
e ^ —U/2), especially for | J| > Ik- see the hollow points 
in Fig. [He)- (d). 

Variation of the conductance with the anisotropy con- 
stant D depends on the sign of exchange parameter J. 
The conductance curves for J > and J < differ signif- 
icantly only for D <^ \J\/2, while when D exceeds | J|/2, 
the difference becomes insignificant. Moreover, it should 
be noted that in the parallel configuration the conduc- 
tance for J > decreases with growing D, whereas for 
J < 0, one observes the opposite tendency, see Fig. [7] 
Such an overlap of the conductance curves for ferromag- 
netic and antiferromagnetic J-coupling in the limit of 
large anisotropy constant D can be explained in a simi- 
lar way as above. One has to take into account that in the 
present situation the difference in energy gaps between 
ground states for J > and J < diminishes, and so do 
the energies of these states as D increases. 

C. Tunnel magnetoresistance (TMR) 

A quantity that describes difference between trans- 
port properties in the parallel (P) and antiparallel (AP) 
magnetic configurations is the tunnel magnetoresistance 
(TMR), defined here asZS 

TMR = ,g . (19) 
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Figure 7. (Color online) Total linear conductance g = g„ 
in the parallel (P - filled points) and antiparallel (AP - hollow 
points) magnetic configuration shown as a function of the 
uniaxial anisotropy constant D for indicated values of the 
J-coupling and two representative OL energies: (a)-(b) e = 
— U/2 and (c)-(d) e = — f//4. Left panel corresponds to the 
ferromagnetic J-coupling (J > 0), while right panel to the 
antiferromagnetic one (J < 0). Other parameters as in Fig. [2] 

Using the conductance data analyzed in the previous 
subsection, we consider now how TMR depends on the 
key parameters of the model, i.e. on the J-coupling, 
Figs.[8^a)-(b) and the magnetic anisotropy D, Figs.[8l^c)- 
(d). 

Let us first discuss the behavior of TMR as a func- 
tion of the exchange coupling constant J, Figs. [Ha)-(b). 
Since for | J| ^ Tk the conductance in the Kondo regime 
is generally larger in the antiparallel magnetic configu- 
ration than in the parallel one, regardless of the type of 
the J-coupling, the corresponding TMR is negative in 
almost the entire Coulomb blockade region. The only ex- 
ception occurs around the particle-hole symmetry point 
e — ~U/2i^ This behavior follows from the suppression 
of the Kondo effect in the parallel configuration due to 
the exchange field, except for the particle-hole symmetric 
point. However, as | J| becomes larger than Tk, the Kondo 
peak becomes suppressed also in the antiparallel configu- 
ration and positive TMR may be observed in the blockade 
regime as well. Moreover, suppression of the conductance 
g^^ for the antiparallel alignment is more evident in the 
case of antiferromagnetic (J < 0) coupling, Fig.[8l[b), and 
the corresponding TMR is therefore significantly larger 
than for the ferromagnetic coupling (J > 0). Another 
observation for J < is that when | J| > Tk, two dis- 
tinctive local maxima develop in the Coulomb blockade 
regime. Their positions depend on J and are symmetri- 
cal with respect to e = —U/2. In addition, for J < 0, the 
TMR considerably surpasses the relevant Julliere's value. 
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Figure 8. (Color online) Tunnel magnetoresistance (TMR) 
presented as a function of the OL energy e and the J-coupling 
(a)-(b) and the magnetic uniaxial anisotropy D (c)-(d). Spe- 
cific data points and parameters used in (a)-(b) correspond to 
those in Fig. [5] and, analogously, in (c)-(d) to Fig. (6) Dashed 
lines are introduced in order to facilitate comparison between 
TMR scales of adjacent plots for J > and J < 0. 



2P^/ (1 — P^)^ which for the present parameters yields 
2/3. Nonetheless, in the regions corresponding to empty 
or doubly occupied OL, one always observes > g^^ 
with TMR approaching the JuUiere's valucj^i 

Consider now variation of TMR with the magnetic 
anisotropy _D, Figs. [5Jc)-(d). First of all, when the OL 
is either empty or occupied by two electrons, TMR re- 
mains insensitive to any changes in the anisotropy con- 
stant D. The same cannot be said about the region of 
£ corresponding to single occupation of the OL, where 
large variation of TMR appears especially for J < 0, 
see Fig. [Hd). From Fig. Efd) follows that the smaller 
the magnetic anisotropy, the larger TMR. Positions of 
the two local maxima in TMR, however, are now rather 
independent of D. Furthermore, for D < |J|/2, TMR 
stays positive. On the contrary, for the ferromagnetic J- 
couphng (J > 0), Fig. [5]^c), the TMR is negative in the 
considered range of e and is only weakly modified upon 
changing D. 



D. The restoring effect of magnetic field 

In the light of the foregoing discussion, we know that 
the Kondo effect is suppressed by exchange field gener- 
ated by ferromagnetic electrodes as well as by the ex- 
change coupling of the OL to magnetic core. Very re- 
cently, it was shown both experimentally and theoreti- 
cally that one can compensate for the exchange- induced 
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Figure 9. (Color online) Restoration of the Kondo peak in 
the total spectral function a(cj) = 1x^2^^ <^^'^i'-^) {^^'^)^ ^-nd 
total conductance g = J]]^ (b,d) by an external magnetic 
field Bz, shown for indicated values of the antiferromagnetic 
J-coupling (J < 0) and for e = —U/2. Left panel corresponds 
to the antiparallel (AP) magnetic configuration, while right 
panel to the parallel (P) one. The other parameters as in 
Fig. [3 



splitting of the orbital level by fine-tuning an external 
magnetic field, restoring thus the universal features of the 
Kondo effect Now we will thus consider the possibility 
of restoring the Kondo effect by applying a compensat- 
ing external magnetic field, B^, oriented along the MQD's 
easy axis, Fig. [51 The interesting observation worth men- 
tioning is that the restoration can take place not only for 
the parallel magnetic configuration. Figs, [nilc)-(d), but 
also for the antiparallel one. Figs. [9l[a)-(b). Furthermore, 
we would like to emphasize that the considered effect oc- 
curs only for the antiferromagnetic (J < 0) coupling be- 
tween an electron in the OL and the MQD's magnetic 
core. From Figs, [ni^c)-(d) follows that for the parallel 
magnetic configuration the full unitary Kondo resonance 
(with the maximum value of conductance g = 1) is re- 
stored, whereas in the antiparallel configuration the peak 
of height 1 — is retrieved owing to the magnetic field. 

From the experimental point of view, it would be use- 
ful to know, at least roughly, how the magnitude of the 
compensating field Be depends on the system's param- 
eters. For this purpose, it is essential to know physical 
mechanism responsible for the restoration of the Kondo 
efi'ect. Thus, in order to gain some deeper understand- 
ing of the problem once again we employ the model of 
a free-standing MQD. It appears that the restoration of 
the Kondo effect becomes possible always when the mag- 
netic field Bz brings the ground state of the spin multiplet 
5t = 5* — 1/2 into resonance with the ground state of the 
multiplet St — S+1/2. Depending on the direction of the 
field, the resonance takes place either between the states 
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) or 1+ |)~ ^ 1+ !)■ Moreover, the possi- 
bility of such degeneracy is straightforwardly granted by 
the fact that these states are characterized by different 
numbers 5* (different zth component of the MQD's to- 
tal spin). Accordingly, the Zeeman contributions to these 
states are different. Because the J-coupling is of the an- 
tiferromagnetic type, the state with greater | S** | has then 
larger energy in the absence of magnetic field {B^ = 0). 
Thus, ii Bz > 0, the energy of the state | — |) decreases 
faster with increasing the field than the energy of | — |)~, 
so the degeneracy of these states occurs at a certain value 
of magnetic field, B ^ Be, which can be then determined 
from the condition A^^m^^^-j ^ g, Fig.[10i;a). 

Taking into account the resonance condition intro- 
duced in the previous paragraph, we find the general ex- 
pression describing dependence of the compensating field 
Be on the J-coupling and magnetic anisotropy D in the 
form 



e^j.,S^oin,i.,n. ,20, 

The comparison between the analytical solution and nu- 
merically derived values of the compensating field B^. for 
different magnetic configurations of the system is pre- 
sented in Fig. [TOlb). It is clearly visible that the nu- 
merical values of B^ generally follow the trend of the 
approximate analytical curve, Eq. (|20p . and substantial 
discrepancies arise only for small values of |J|, as one 
might expect. The slightly smaller value of the compen- 
sating field in the situation of a MQD attached to mag- 
netic electrodes can be attributed to renormalization of 
MQD's energy levels due to the strong OL-electrodes cou- 
pling, which results in diminishing energy gaps between 
the states participating in formation of the Kondo effect. 
It should be here emphasized that this difference can be 
seen only due to specific normalization of the compen- 
sating field Be with respect to \J\. Otherwise, when pre- 
sented in the logarithmic scale, the curves correspond- 
ing to numerical and analytical solutions follow the same 
trend, see the inset in Fig. [TUTb). 

Finally, since the form of Eq. (PD|) is exactly the same 
as for A^^'^, Eq. (|B2p . the asymptotic values of Be are 
immediately obtained as: Be ~ (25*4-1)1 J|/2 for | J| > D, 
and Be ~ S\J\ for |J| ^ D. Additionally, it might be 
helpful to know how the change of magnetic anisotropy 
D infiuences the analytical solution Be{\J\), bold line in 
Fig. [TUr bV For this purpose, we calculate the inflexion 
point position (IPP) of the compensating magnetic field 
curve Be{\J\), 



IPP 



2(2g- 1) 
(25+1)2 



l-{2S)^/^]\l + i2S)^^^]D. (21) 



It turns out that IPP depends linearly on D, with the 
proportionality constant being a complex function of the 
MQD's core spin number S. Consequently, modification 
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Figure 10. (Color online) (a) Energies of ground states | — §) 
and |— §)~ shown as functions of an external magnetic field 
for indicated values of the exchange coupling parameter J and 
for D/U ^ 1.7-10-* {D/Tk = 0.1) and e = -[7/2. Dashed line 
corresponds to the analytical solution given by Eq. (|20|) . (b) 
Dependence of the compensating field Be on the value of the 
J-coupIing. Points represent numerical results obtained for 
the antiparallel (AP) and parallel (P) magnetic configuration, 
while the bold line delineates the analytical solution (a.s.), 
Eq. (|20|) . IP indicates the inflexion point of the curve Bc(| J|). 
Inset: the compensating field Be normalized to Tk, shown in 
the logarithmic scale. Other parameters as in Fig. (2] 



of the magnetic anisotropy constant D does not affect the 
general shape of the compensating field curve Be{\J\), 
but it only leads to shifting of the curve's inflexion point. 



IV. SUMMARY AND CONCLUSIONS 

In this paper we have investigated transport properties 
in the Kondo regime of a class of systems exhibiting uni- 
axial magnetic anisotropy. The model assumed includes 
one orbital level through which electrons can tunnel and 
which is additionally exchange-coupled to a magnetic mo- 
ment. The model can describe a single-level quantum dot 
in which electrons are exchange-coupled to an embedded 
magnetic impurity. It can also be used to describe trans- 
port through magnetic atoms and molecules. 

Using the numerical renormalization group method we 
have calculated spectral density of the OL level and linear 
conductance of the system. The key new feature of trans- 
port characteristics is the suppression of the Kondo effect 
by exchange coupling to magnetic core. Independently of 
the sign of the coupling parameter J, suppression takes 
place for both ferromagnetic as well as antiferromagnetic 
coupling. In the limit of small J, | J| -C Tk, we find the 
Kondo resonance characteristic of single-level dots. This 
resonance is suppressed with increasing J but then some 
side resonances appear in the spectral density. It is wor- 
thy of note that the suppression of the Kondo peak in 
the case considered appears gradually with increasing J, 
contrary to the case of a dot exchange-coupled to elec- 
tron reservoir of continuous density of states, where the 
suppression is associated with a quantum phase transi- 
tion^ 

We have also shown that the suppressed Kondo reso- 
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nance in spectral function and in transport characteris- 
tics can be restored by application of an external mag- 
netic field. This restoration is however complete only in 
the case of antiferromagnetic exchange coupling ( J < 0) 
for both the parallel and antiparallel magnetic configu- 
rations. 
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Appendix A: Calculation of the spectral function 

During the course of NRG iterative calculation, the 
spectral function Acr{uj) is obtained as a set of the Dirac 
delta functions d{uj — LUn) at frequencies w„, which next 
have to be broaden in order to acquire a continuous spec- 
trum. Because of logarithmic discretization of conduction 
band, it is convenient to collect the delta peaks in loga- 
rithmic bins and broaden them using a logarithmic Gaus- 
sian distribution, with the broadening parameter typi- 
cally given by log(A)i^i^i2^ It turns out, however, that 
due to logarithmic discretization of the band and trunca- 
tion during the NRG run, the broadened spectral func- 
tion may exhibit some artifacts, such as an oscillatory 
behavior for energies smaller than the Kondo tempera- 
ture. One of the tricks to overcome these problems is the 
so-called self-energy trickS^i^ 

The essential idea of the self-energy trick relies on find- 
ing the spectral function by constructing the full self- 
energy Ecr(w) of the system)^ The self-energy can be 
expressed as a ratio of two spectral functions. Although 
each spectral function displays similar problems related 
with discretization, by calculating their ratio one obtains 
a smooth function where most of the artifacts are sup- 
pressed. Having found the self-energy, one can then calcu- 
late the retarded Green's function ((Cg.|cJ.))[, of the OL, 
Eq. (fTT|) . using the equation of motion method^ 



((c.l4))L = 



1 



(Al) 



where Eo- = e -f ^Bz{5a\ — ^ai)- The total self-energy 
consists of three terms, 

= E^'-(^) + + A:,(c.), (A2) 

which represent contributions stemming from: the 
Coulomb interaction, E^''(w), the exchange interaction 



Finally, the improved spectral function of the orbital level 
is given by 



1 



ImE^(^) 



^ [uj-ea-Rei:i{uj)\ +[lmS',H]' 
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Appendix B: Energy gaps 

In the case of the ferromagnetic (FM) J-coupling (J > 
0) , exact analytical expressions for the energy gaps A™ 
and A™ can be derived as the energy difference between 
appropriate MQD's states: |± |) and |± for Af'^, 
and I ± I) and I ± §)+ for l^l^p. so that 



aFM _ 2Sjfl 2S 1 

^1(2) - . \'J\ + 7) 



D 



16 



The energy gaps for antiferromagnetic (AFM) J-coupling 
( J < 0) can be found in a similar way. The energy differ- 
ence between the MQD's states |± |)~ and |± |) is then 



given by 



^AFM ^ ^ + 1 I T, 25-1 



-\J\ 



D 



16 



whereas the formula for the gap between the states |±|)" 



and |± |)+ takes the form 



^AFM ^ 2a/I?^^±11!j2 + ^^l^D{D + I J|). 
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